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ABSTRACT 

Direct  numerical  simulation  data  bases  for  compressible  homogeneous  shear  flow  are 
used  to  evaluate  the  performance  of  recently  proposed  Reynolds  stress  closures  for  com¬ 
pressible  turbulence.  Three  independent  pressure-strain  modeb  are  considered  along  with 
a  variety  of  explicit  compressible  corrections  that  account  for  dilatational  dissipation  and 
pressure-dilatation  effects.  The  ability  of  the  models  to  predict  both  time  evolving  fields  and 
equilibrium  states  is  systematically  tested.  Consutent  with  earlier  studies,  it  is  fotmd  that 
the  addition  of  simple  dilatational  modeb  allows  for  the  prediction  of  the  reduced  growth 
rate  of  turbulent  kinetic  energy  in  compressible  homogeneous  shear  flow.  However,  a  closer 
examination  of  the  equilibrium  structural  parameters  uncovers  a  major  problem.  None  of  the 
modeb  are  able  to  predict  the  dramatic  increase  in  the  normal  Reynolds  stress  anisotropies 
or  the  significant  decrease  in  the  Reynolds  shear  stress  anisotropy  that  arise  from  compress¬ 
ible  effects.  The  physical  origin  of  this  deficiency  is  attributed  to  the  neglect  of  compressible 
terms  in  the  modeling  of  the  deviatoric  part  of  the  pressure-strain  correlation. 


*This  leseaich  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NASl-19480  while  the  author  was  in  rendence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681-0001. 


1.  INTRODUCTION 


The  need  to  predict  supersonic  and  hypersonic  turbulent  flows  of  aerodynamic  importance 
has,  in  recent  years,  led  to  new  reseuch  initiatives  in  compressible  turbulence  modeling.  A 
major  stumbling  block  in  the  development  of  improved  compressible  turbulence  models  is  the 
lack  of  detailed  experimental  data  for  the  compressible  turbulence  statistics  in  basic  high¬ 
speed  compressible  flows.  Elxperimental  limitations  currently  make  it  infeasible  to  obtain 
detailed  measurements  of  any  turbulence  statistics  beyond  the  mean  velocity  and  Reynolds 
shear  stress  in  supersonic  turbulent  flows.  This  makes  it  virtually  impossible  to  pinpoint  the 
origin  of  deficient  model  predictions  when  they  arise.  In  Reynolds  averaged  calculations  of 
complex  supersonic  turbulent  flows,  erroneous  predictions  for  the  mean  velocity  field  may 
arise  from  modeling  errors  in  the  Reynolds  stresses  that  can  be  traced  to  a  variety  of  possible 
deficiencies  in  the  treatment  of  compressibility  effects. 

During  the  past  four  years,  direct  numerical  simulation  (DNS)  data  bases  of  some  basic 
supersonic  turbulent  flows  have  become  available  [1, 2],  wherein  the  full  compressible  Navier- 
Stokes,  continuity  and  energy  equations  are  solved  numerically  with  all  turbulent  scales 
resolved.  Most  notable  among  these  are  the  DNS  data  bases  for  compressible  homogeneous 
shear  flow  (see  Blaisdell,  Mansour  and  Reynolds  [3]  and  Sarkar,  Erlebacher  and  Hussaini 
[4,5]).  Homogeneous  shear  flow  is  an  important  building-block  flow  since  it  encapsulates 
some  of  the  important  features  of  an  equilibrium  turbulent  boundary  layer  in  a  simplified 
setting  that  is  unencumbered  by  the  effects  of  turbulent  diffusion  or  wall  blocking.  It  has 
been  an  extremely  useful  test  case  for  the  calibration  and  screening  of  turbulence  models  for 
incompressible  flows  [6-8].  However,  no  comprehensive  test  of  compressible  Reynolds  stress 
models  has  been  made  using  these  new  DNS  data  bases  for  compressible  homogeneous  shear 
flow.  This  establishes  the  motivation  for  the  present  paper.  Prior  to  this  study,  only  very 
limited  comparisons  of  compressible  turbulence  modeb  have  been  made  using  these  DNS 
results  [9,10]. 

In  this  paper,  the  DNS  data  base  of  Blaisdell,  Mansour  and  Reynolds  [3]  for  compressible 
homogeneous  shear  flow  will  be  used  to  critically  evaluate  recently  proposed  compressible 
Reynolds  stress  closures.  Full  second-order  closures  will  be  considered  that  explicitly  ac¬ 
count  for  high-speed  compressible  effects.  Variable  density  extensions  of  three  incompress¬ 
ible  pressure-strain  models  due  to  Launder  and  co-workers  [11,  12]  and  Speziale,  Sarkar  and 
Gatslri  [13]  will  be  tested  along  with  models  for  the  compressible  dissipation  and  pressure- 
dilatation  that  were  recently  proposed  by  Sarkar  et  ai  [14,15]  and  Zeman  [16,17].  The  ability 
of  each  model  to  predict  time-evolving  fields  and  equilibrium  states  such  as  the  Reynolds 
stress  anisotropies  will  be  assessed  in  detail.  An  attempt  will  also  be  made  to  understand 
the  physical  origin  of  deficient  model  predictions.  In  this  regard,  some  important  results  will 
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be  uacoveied  that  were  overlooked  in  previont  stndiei.  These  issues  will  be  discussed  fully  in 
the  sections  to  follow  and  a  recommendation  will  be  made  for  the  development  of  improved 
modeb. 

2.  THEORETICAL  BACKGROUND 

We  will  consider  the  mean  turbulent  flow  of  an  ideal  gas  governed  by  the  compressible 
Navier-Stokes  equations.  The  mass  density  p  is  decomposed  into  standard  ensemble  mean 
and  fluctuating  parts,  respectively,  as  follows: 


p  =  p  +  /• 


(2.1) 


The  velocity  field  v,-  and  temperature  T  are  decomposed  into  standard  ensemble  mean  and 
fluctuating  parts  |pven  by 

tt,  =  ti,  +  ttj,  T  =  T  +  T'  (2.2) 

or  into  mass- weighted  mean  and  fluctuating  parts  pven  by 


u,  =  tt,  -H  uf ,  r  =  T  +  T". 


(2.3) 


For  any  flow  variable  the  quantity  T  =  ^/p  denotes  the  Favre  or  mass- weighted  average 
given  that  an  overbar  represents  a  standard  ensemble  mean. 

The  Favre-averaged  continuity,  Navier-Stokes  and  energy  equations  take  the  form  [18, 
19] 

^  +  (?Si)i  =  0  (2.4) 

Q 

^(^)  +  (?*»i«i)a  =  (2.5) 

-I-  ¥  (2.6) 

in  homogeneous  turbulent  flows  where  the  mean  velocity  gradients  and  all  higher-order  statis¬ 
tics  are  spatially  uniform.  In  Eqs.  (2.5)  and  (2.6), 


r  =  pier 


(2.7) 


*  = 

“■«  =  -  +  “w) 


(2.8) 

(2.9) 


are,  respectively,  the  thermodynamic  pressure,  viscous  dissipation  and  viscous  stress  tensor; 
R  is  the  ideal  gas  constant,  Cv  is  the  specific  heat  at  constant  volume,  p  is  the  dynamic 
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viscodty  of  the  gai,  end  (  )^'  =  d(  )/dst  denotes  a  spatial  gradient.  Since  the  tuibnlent  mass 
flnx  vanishes  in  onnpressible  homogeneous  turbulence,  it  fidlows  that  [19] 


¥  =  VijUij  +  (2.10) 

at  high  turbulence  Reynolds  numbers  where  e  =  is  the  turbulent  dissipation  rate. 

Following  the  recent  work  of  Sarkar  et  al  [14]  and  Zeman  [16],  the  turbulent  dissipation 
rate  can  be  decomposed  into  solenoidal  and  dilatational  parts  as  follows: 

e  =  e.  +  Cc  (2.11) 

where,  for  homogeneous  turbulence, 

e-  =  ec  =  |/*(tt|,,)V?  (212) 

are,  respectively,  the  solenoidal  and  compressible  (or  dilatational)  parts  of  the  turbulent 
dissipation  rate  pven  that  is  the  fluctuating  vortidty.  Here,  e.  represents  the  turbulent 
dissipation  arising  from  the  traditional  ene^y  cascade  which  is  solenoidal  (i.e.,  vortical) 
in  character;  £«  represents  the  turbulent  distipation  arising  from  purely  compressible  or 
dilatational  modes  (««  =  0  for  incompressible  turbulent  flows).  For  homogeneous  turbulence, 
this  decomposition  is  unique. 

Since  the  mean  pressure 

f  =  pBT,  (2.13) 

it  is  clear  from  (2.4)-(2.10)  that  closure  of  the  mean  flow  equations  is  achieved  in  compressible 
homogeneous  turbulence,  at  high  Reynolds  numbers,  once  models  for  the  pressure-dilatation, 
compressible  dissipation  and  solenoidal  distipation  are  provided.  The  determination  of  these 
correlations  also  requires  information  on  the  turbulence  intensity  level  and,  hence,  on  the 
Reynolds  stress  tensor.  The  Favre-averaged  Reynolds  stress  tensor  Tij  =  u"u"  is  a  solution 
of  the  transport  equation  [18, 19] 

2  2 _ 

pTij  =  -pTikUjjk  -  +  Ho  -  -peSij  -h  (2.14) 

which  is  exact  for  compressible  homogeneous  turbulence  pven  that 

n<i  =  + 

- 

is  the  difference  between  the  deviatoric  parts  of  the  pressure-strain  correlation  and  distipation 
rate  tensor.  Eq.  (2.14)  contains  the  pressure-dilatation  correlation  p't^  as  well  as  the 


compiesaible  mnd  lolenddal  parts  of  the  dissipation  rate  tensor  dnce  e  =  e*  +  Cc.  Hence, 
a  full  Reynolds  stress  closure  is  achiered  in  compressible  turbulence  if  models  are  provided 
for: 

(i)  The  difference  between  the  deviatoric  parts  of  the  pressure-strain  correlation  and  dis¬ 
sipation  rate  tensor,  Hij. 

(ii)  The  pressure-dilatation  correlation  p'i^. 

(iii)  The  solenoidal  dissipation  e.. 

(iv)  The  compressible  dissipation  Ce. 

In  compressible  homogeneous  shear  flow,  the  mean  velocity  gradient  tensor  is  given  by 

Uij  =  S6iiSj2 

where  S  is  the  shear  rate  which  is  constant.  The  mean  density  is  constant  whereas  the  mean 
temperature  is  a  function  of  time  alone,  i.e., 

p  =  constant  (2.17) 

f  =  f  (t).  (2.18) 

E^s.  (2.16)  and  (2.17),  which  are  the  same  as  their  incompressible  counterparts,  identically 
satisfy  the  mean  continuity  and  mean  momentum  equations  (2.4)  and  (2.5).  Assuming  that 
the  mean  spedflc  heat  is  constant,  the  Reynolds-averaged  energy  equation  (2.6)  simplifies 
to  the  form 

jCvT  =  +  pe  (2.19) 

for  a  compressible  homogeneous  shear  flow  at  high  Reynolds  numbers. 

The  substitution  of  (2.16)  into  the  contraction  of  (2.14)  yields  the  turbulent  kinetic  energy 
equation 

pk  =  fP-¥f^ji-pe  (2.20) 

where  K  =  is  Favre-averaged  turbulent  kinetic  energy  and  V  =  —tx^S  is  the 

turbulence  production.  Equation  (2.20)  can  be  combined  with  (2.19)  to  yield  a  transport 
equation  for  the  turbulence  Mach  number.  This  transport  equation  takes  the  form 

M,  =  ^7>  +  ^[l  +  i7(7-l)«?](pX;-?e)  (2-21) 

where  7  =  is  the  ratio  of  specific  heats  and  Mt  =  yj2KI'iBT  is  the  turbulent  Mach 

number.  Recently  proposed  models  for  the  four  turbulence  correlations  that  are  needed  for 
closure  -  namely,  e„  and  Hij  -  will  be  discussed  in  the  next  section. 
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S.  THE  TURBULENCE  MODELS  TO  BE  TESTED 


In  almost  aU  compressible  second-order  closures,  the  deriatoiic  part  of  the 

pressure-strain  correlation  is  modeled  as  a  variable  density  extension  of  its  incompressible 
counterpart.  FSirthermore,  as  in  the  majority  of  incompressible  turbulence  models,  the  devi- 
atoric  part  of  the  dissipation  rate  tensor  is  neglected  by  invoking  the  Kolmogorov  assumption 
of  local  isotropy  [18, 19].  This  leads  to  a  hierarchy  of  models  that  are  of  the  general  form: 

+  pKMiju{b){UkJ  —  ^Um,mSu)  (3.1) 


where 


hij  = 


« 


is  the  Reynolds  stress  anisotropy  tensor.  In  (3.1),  Aij  and  Miju  are  identicai  to  their  in¬ 
compressible  forms;  compressibility  effects  are  only  accounted  for  through  changes  in  the 
mean  density.  Since  the  mean  density  Jl  is  constant  in  homogeneous  shear  flow  -  and  the 
mean  dilatation  u,-,,-  is  aero  -  it  follows  that  the  hierarchy  of  pressure-strain  models  (3.1) 
is  completely  identical  to  its  incompressible  counterpart.  The  consequences  of  this  will  be 
discussed  later. 


In  this  paper,  variable  density  extensions  of  three  incompressible  pressure-strain  modeb 
will  be  considered:  the  Launder,  Reece  and  Rodi  [11]  model,  the  Fu,  Launder  and  Tselepi- 
dakis  [12]  model  and  the  Speaiale,  Sarkar  and  Gatsld  [13]  model.  The  first  model  is  chosen 
since  it  is  the  most  widely  used  pressure-strain  model;  the  last  two  models  are  chosen  since 
they  have  recently  been  shown  to  perform  the  best  among  a  variety  of  existing  models  for  in¬ 
compressible  homogeneous  shear  flow  [7,8].  The  detailed  form  of  these  modeb  are  provided 
below. 


Launder,  Reece  &  Rodi  Model 

Ilij  =  —Cipe,bij  +  -pK  ^Sij  —  -SkkSi^  +  C2pK 

-  2  -  \  _ 

+bjkSik  —  •^buSuSijj  +  C3pK{bik<tfjk  +  bjifi/ik) 

where  ^ 

Sij  =  2(®*iJ  *it»)*  ~  2^^'^  ~ 

Cl  =  3.0,  C2  =  1.75,  C3  =  1.31 


(3.3) 

(3.4) 

(3.5) 
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At,  JjMmder  &  Ttekpidakia  Modd 

n<i  =  —C\pe,hii  +  CiPe,  -  ^bybuSij^ 

(sy  -  + 1.2?* 

+bitSik  —  ^htiSuSi^  + 

+6jfcw<fc)  +  ^pK(bikbuSjt  +  bjkbuSa 

—2bik^ubij  —  3bul>ubij) 

4_  -  .14 

+^P^(bitbuuiji  +  bjkbuufii)  —  —pKl8II(bijtu>jk 

-hbjkibik)]  +  12(biitbuunmbm3  + 

(3.«) 

where 

Cl  =  2  -  120//A*/’  -  Ca  =  144//A»/* 

(3.7) 

(3.8) 

F=  1  +  K//  +  2T//I 

(3.9) 

Speziak,  Sarkar  &  Gatshi  Mode/ 

Dy  =  —{Cipe,  4-  CipP)bij  +  C2pe,  ^^bkj  —  ^bubySi^ 

+(173  -  C;lll'')pK  (Sii  -  \Sufii) 

+C4?*  ^6135^  +  h^Sik  —  zbuSubi^ 

+CspK(biifijk  +  bjkufik) 

(3.10) 

where 

Cl  ~  3.4,  Cl  =  1.80,  Ca  —  4.2 

(3.11) 

Ca  =  1,  C;  =  1.30,  <74  =  1.25 

5 

(3.1S) 

Ci  =  0.40,  Ih  =  bijbij^  V  =  —TiiUij. 

(3.13) 

These  models  will  hereafter  be  referred  to  as  the  LRR,  FLT  and  SSG  models,  respectively. 
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Since  the  mean  dilatation  n,.,-  ii  aeto  in  homogeneoui  ihear  flow,  virtually  all  eziating 
model  tiamport  equationa  for  the  aolenmdal  diaaipation  are  of  the  form  [19  -  21] 

i.  =  C.t^r-Ca^  (3.14) 

in  compreaaible  homogeneoua  ahear  flow.  Thia  equation  ia  idetUical  to  ita  incompreaaible 
counterpart.  In  (3.14),  V  =  —rx^S  ia  the  turbulence  production  and  Cei  and  Cg2  ue  con- 
atanta  (in  the  SSG  model,  C,x  —  1.44  and  Cg2  =  1.83  whereaa  in  the  LRR  and  FLT  modela, 
Cgi  =  1.44  and  Cca  =  1.90).  Two  recent  modela  for  the  compressible  dissipation  will  be  con¬ 
sidered  that  were  proposed  by  Sarkar  et  al  [14]  and  Zeman  [16].  These  models  are  algebraic 
and  of  the  general  form 

ec  =  (3.15) 

where  Mt  is  the  turbulence  Mach  number  defined  earlier.  In  the  Sarkar  et  al  [14]  model 

Ce  =  otiMfe,  (3.16) 


which  is  obtained  from  an  asymptotic  analysis  that  is  formally  valid  for  small  turbulent  Mach 
numbers.  The  constant  ai  in  (3.16)  was  determined  to  be  approximately  0.5  based  on  direct 
numerical  simulations  of  homogeneous  turbulence  [15].  On  the  other  hand,  the  compressible 
dissipation  rate  model  of  Zeman  [16]  is  based  on  an  analysis  that  incorporates  the  effects  of 
eddy  shocklets.  For  homogeneous  shear  flow,  this  model  takes  the  form  [20] 

Cc  =  {l  -  exp(-[(J»f,  -  0.25)/0.8]*)}  c,  (3.17) 

for  Mt  >  0.25;  Cc  =  0  for  <  0.25. 

The  pressure-dilatation  model  that  we  will  primarily  consider  ia  that  due  to  Sarkar  [15]. 
This  model  is  algebraic  and  takes  the  following  simple  form  in  homogeneous  shear  flow: 

f'u'ii  =  -a^fPMt  +  a3pe,M^  (3.18) 


where  02  and  as  are  constants  that  take  on  the  values  of  0.15  and  0.2,  respectively.  Some 
limited  comparisons  will  also  be  made  with  the  pressure-dilatation  model  of  Zeman  [17] 
which  is  (pven  by 


(3.19) 


where 


(3.20) 
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and 


Of  course,  in  (3.19),  p  s  'pRT  and  u,,,  =  0  for  homogeneous  shear  flow.  In  the  next  section, 
a  comparison  of  computed  results  obtained  from  these  modds  for  compressible  homogeneous 
shear  flow  will  be  made. 

4.  DISCUSSION  OF  RESULTS 

The  transport  equations  (2.6),  (2.14)  and  (3.14)  -  incorporating  the  models  discussed  in 
Section  3  -  were  solved  numerically  for  compresdble  homogeneous  shear  flow  using  a  fourth- 
order  accurate  Runge-Kutta  numerical  integration  scheme.  Comparisons  will  be  made  with 
the  direct  numerical  simulations  (DNS)  of  compresnble  homogeneous  shear  flow  by  Blaisdell 
et  al  [3].  Both  time-evolving  fields  and  equilibrium  states  will  be  compared  with  DNS  results. 
Comparisons  will  be  made  with  run  SHA  192  of  Blaisdell  et  al  [3]  since  it  is  the  best  resolved 
run  that  is  of  long  duration  (until  St  »  24).  It  will  be  shown  that,  consistent  with  the  DNS, 
the  models  predict  that  the  turbulent  kinetic  miergy,  turbulent  dissipation  rate  and  mean 
temperature  grow  exponentially  with  time.  All  structural  parameters  such  as  the  anisotropy 
tensor  bij  and  turbulent  Mach  number  Af|  achieve  equilibrium  values  that  are  independent 
of  the  initial  conditions. 

First,  we  will  consider  model  predictions  for  the  case  where  there  are  no  explicit  dilata- 
tional  terms.  This  is  done  for  one  major  reason:  until  very  recently,  the  vast  majority  of 
Reynolds  stress  calculations  of  compressible  turbulent  shear  flows  were  conducted  with  vari¬ 
able  density  extensions  of  incompressible  models  where  explicit  turbulent  dilatational  terms 
were  not  included.  Consequently,  a  benchmark  is  established  for  assessing  the  performance 
of  the  recently  proposed  dilatational  models.  In  Figure  1,  the  time  evolution  of  the  turbulent 
kinetic  energy  predicted  by  variable  density  extensions  of  the  LRR,  FLT  and  SSG  models  - 
with  compressible  dissipation  and  pressure-dilatation  effects  neglected  -  are  compared  with 
DNS  results  of  Blaisdell  et  al  [3]  (run  SHA  192).  Here,  AT*  =  KfKo  and  t*  =  St  are  the 
dimensionless  turbulent  kinetic  energy  and  the  dimennonless  time,  respectively,  where  a 
subscript  0  denotes  the  initial  value.  fVom  these  results,  it  is  clear  that  all  of  the  models 
drastically  overpredict  the  growth  rate  of  the  turbulent  kinetic  energy  -  a  deficiency  that 
arises  from  the  complete  ne^ect  of  turbulent  dilatational  effects  as  discussed  by  previous 
authors  [14,  15,  19].  In  Figures  2-3,  the  model  predictions  for  the  time  evolution  of  the 
turbulent  dissipation  rate  (e*  =  e/eo)  and  turbulent  Mach  number  (Jlf|)  are  displayed.  Here, 
as  with  the  results  shown  in  Figure  1  for  the  turbulent  kinetic  energy,  the  calculations  are 
made  for  >  2  in  order  to  avoid  the  unphysical  early  transient  of  the  DNS  which  lasts  for 
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appmodmatelj  one  eddy  inmover  time  due  to  the  fact  that  the  flow  is  initially  seeded  with 
an  artifidal  random  Gaussian  fleld.  It  is  obvious  from  these  results  that  the  models  substan¬ 
tially  overpredict  the  turbulent  dissipation  rate  and  turbulent  Mach  number.  In  regard  to 
the  latter,  highly  unphysical  equilibrium  turbulence  Mach  numbers  greater  than  1.5  are  pre¬ 
dicted.  This  is  true  for  both  the  FLT  and  SSG  models  which  yield  excellent  predictions  for 
incompressible  homogeneous  shear  flow.  Consequently,  it  is  clear  that  variahlt  density  exten¬ 
sions  of  incompressibk  Reynolds  stress  models,  with  turhtdeni  dilatational  effects  neglected, 
cannot  properly  describe  compressible  homogeneous  shear  flow  -  a  conclusion  consistent  with 
previous  findings.  The  compressible  dissipation  and  pressure-dilatation  correlation  pve  rise 
to  a  significant  reduction  in  the  growth  rate  of  the  turbulent  kinetic  energy  for  compressible 
homogeneous  shear  flow. 

In  Figure  4,  computed  results  for  the  turbulent  kinetic  energy  obtained  from  the  LRR, 
FLT  and  SSG  models  -  with  the  compressible  dissipation  and  pressure-dilatation  models 
of  Sarkai  et  al  [14,  15]  -  are  compared  with  DNS  results.  With  the  addition  of  these 
turbulent  dilatational  terms,  the  FLT  and  SSG  models  are  now  able  to  properly  predict 
the  reduced  growth  rate  of  the  turbulent  kinetic  energy  that  arises  from  compressibility 
effects  in  homogeneous  shear  flow.  The  LRR  model  overpredicts  the  growth  rate  of  the 
turbulent  kinetic  energy  by  an  amount  comparable  to  that  which  has  been  documented  for 
incompressible  homogeneous  shear  flow  (Speziale,  Gatski  and  Mac  Giolla  Mhuiris  [7]  showed 
that  the  LRR  model  overpredicts  the  growth  rate  by  about  25%  for  the  incompressible 
case).  The  computed  time  evolutions  of  the  turbrdent  dissipation  rate  are  compared  with 
DNS  results  in  Figure  5.  Although  the  dissipation  rate  is  overpredicted,  the  inclusion  of  the 
dilatational  modek  of  Sarkar  leads  to  a  substantial  improvement  over  the  results  shown  in 
Figure  2.  The  computed  model  predictions  for  the  time  evolution  of  the  turbulence  Mach 
number  are  compared  with  DNS  results  in  Figure  6.  With  the  addition  of  the  dilatational 
models  of  Sarkar,  the  FLT  and  SSG  models  now  yield  remarkably  good  predictions  for  the 
turbulence  Mach  number  as  illustrated  by  these  results. 

Computed  results  for  the  LRR,  FLT  and  SSG  models  with  the  addition  of  the  compress¬ 
ible  dissipation  and  pressure-dilatation  models  of  Zeman  [16,  17]  will  now  be  considered. 
Model  predictions  for  the  time  evolution  of  the  turbulent  kinetic  energy,  turbident  dissipa¬ 
tion  rate  and  turbident  Mach  number  are  compared  with  DNS  results  in  Figures  7-9.  Most 
notably,  the  predictions  of  the  FLT  and  SSG  models  for  the  turbulent  kinetic  energy  and 
turbulent  Mach  number  are  comparably  good  to  those  obtained  using  the  dilatational  mod¬ 
els  of  Sarkar.  Alternative  tests  of  the  iSeman  models  were  recently  reported  for  compressible 
homogeneous  shear  flow  that  are  comparably  favorable  [22]. 

In  Figures  10(a)- 10(c),  computed  results  for  the  independent  non-vanishing  components 


9 


of  the  Reynolds  stress  anisotropy  tensor  obtained  from  the  LRR,  FLT  and  SSG  models  - 
using  the  dilatational  modds  of  Sarkar  -  are  compared  with  DNS  results.  It  is  clear  from 
these  results  that  all  of  the  models  drastically  underpredict  the  magnitude  of  the  normal 
Reynolds  stress  anisotropies  in  compressible  homogeneous  shear  flow;  on  the  other  hand,  the 
modek  overpredict  the  magnitude  of  the  Reynolds  shear  stress  anisotropy.  The  DNS  results 
of  Blaisdell  et  ai  [3]  indicate  that  the  normal  Reynolds  stress  anisotropies  in  compressible 
homogeneous  shear  flow  are  nearly  twice  as  large  as  their  incompressible  counterparts  [23]; 
the  magnitude  of  the  Reynolds  shear  stress  anisotropy  is  approximately  25%  lower  for  com¬ 
pressible  homogeneous  shear  flow.  Even  with  the  addition  of  what  appear  to  be  reasonably 
sound  models  for  the  compressible  dissipation  and  pressure  dilatation,  variable  density  ex¬ 
tensions  of  the  LRR,  FLT  and  SSG  models  are  unable  to  predict  the  dramatic  changes  in 
the  Reynolds  stress  anisotropies  that  arise  from  compressible  effects.  This  can  be  seen  more 
clearly  from  Tables  1-3  where  the  model  predictions  for  the  equilibrium  Rejmolds  stress 
anisotropies  are  compared  vrith  DNS  results  for  incompressible  homogeneous  shear  flow  and 
for  compressible  homogeneous  shear  flow,  wherein  both  the  dilatational  models  of  Sarkar  and 
Zeman  are  implemented.  It  is  obvious  from  these  results  that  the  individual  model  predic¬ 
tions  are  in  reasonably  dose  range  for  both  compressible  and  incompressible  homogeneous 
shear  flow  whereas  the  DNS  results  are  drastically  different. 

5.  CONCLUSION 

A  systematic  evaluation  of  recently  proposed  Reynolds  stress  turbulence  dosures  for 
high-speed  compressible  flows  has  been  conducted  with  the  use  of  the  DNS  data  base  of 
Blaisdell  et  al  [3]  for  compressible  homogeneous  shear  flow.  The  recently  devdoped  dilata¬ 
tional  models  of  Sarkar  et  al  [14,  15)  and  Zeman  [16,  17]  were  tested  in  conjunction  vdth 
variable  density  extensions  of  two  of  the  newest  modek  for  the  pressure-strain  corrdation: 
the  Fu,  Launder  and  Tsdepidakis  (FLT)  modd  and  the  Speziale,  Sarkar  and  Gatski  (SSG) 
modd.  Consistent  with  the  findings  of  earlier  studies,  it  was  found  that  when  dilatational 
effects  arising  from  the  compressible  dissipation  and  the  pressure-dilatation  corrdation  are 
neglected,  variable  density  extensions  of  existing  Reynolds  stress  turbulence  dosures  yidd 
poor  predictions  for  compressible  homogeneous  shear  flow.  DNS  results  indicate  that  com¬ 
pressibility  effects  lead  to  a  substantial  reduction  in  the  growth  rate  of  the  turbulent  kinetic 
energy  -  a  feature  that  cannot  be  predicted  by  Reynolds  stress  turbulence  dosures  wherein 
turbulent  dilatational  terms  are  neglected.  With  the  addition  of  the  compressible  dissipation 
and  pressure-dilatation  modek  of  Sarkar  and  Zeman,  the  newest  Reynolds  stress  turbulence 
dosures  are  able  to  accuratdy  predict  the  reduction  in  the  growth  rate  of  the  turbulent 
kinetic  energy.  Accurate  results  are  also  obtained  for  the  turbulent  Mach  number  which  is 
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overpzedicted  by  100%  when  toibnlent  dilntational  terms  are  neglected.  There  appears  to 
be  little  doubt  these  dilatational  models  represent  significant  progress  in  the  modeling  of 
compressible  homogeneous  turbulence. 

Although  the  inclusion  of  the  newest  compresnble  dissipation  and  pressure-dilatation 
models  in  Reynolds  stress  turbulence  closures  leads  to  substantially  improved  predictions  for 
the  turbulent  kinetic  energy  and  Mach  number,  some  substantial  deficiencies  still  remain  for 
the  proper  description  of  compressible  hom<^eneous  shear  fiow.  DNS  results  indicate  that 
compressibility  effects  in  homogeneous  shear  flow  lead  to  a  substantial  modification  of  the 
equilibrium  Reynolds  stress  anisotropies  wherein  the  magnitude  of  the  normal  components 
are  nearly  doubled  and  the  magnitude  of  the  shear  component  is  reduced  by  approximately 
25%.  Even  when  compressible  dissipation  and  pressure-dilatation  models  are  added,  the 
Reynolds  stress  turbulence  closures  considered  in  this  study  are  still  unable  to  predict  this 
effect.  This  is  due  to  deficiencies  in  the  modeling  of  the  deviatoric  part  of  the  pressure- 
strain  correlation  which  controls  the  level  of  Reynolds  stress  anisotropy.  The  pressure- 
strain  models  considered  herein  -  as  well  as  those  used  in  virtually  all  previous  studies 
of  compressible  turbulence  -  do  not  account  for  explicit  compressible  effects.  This  is  a 
reasonable  apprommation  for  compressible  flows  where  the  turbulent  Mach  number  Af(  < 
0.3:  a  restriction  that  allows  for  the  Morkovin  [24]  hypothesis  to  be  invoked.  However,  in 
compressible  homogeneous  shear  flow  the  turbulence  Mach  number  achieves  an  equilibrium 
value  of  approximately  0.6  -  a  value  that  is  too  large  to  neglect  explicit  compressible  effects 
in  the  modeling  of  the  pressure-strain  correlation.  The  standard  hierarchy  of  pressure-strain 
models  is  based  on  an  analysis  of  the  incompressible  Poisson  equation  for  the  pressure  that 
does  not  apply  to  high-speed  compressible  flows.  For  these  flows,  the  pressure  is  determined 
from  a  thermodynamic  equation  of  state  unlike  the  pressure  in  incompressible  flows  which 
is  a  Lagrange  multiplier  determined  by  the  solenoidal  constraint  on  the  velocity  field.  If 
high-speed  turbulent  shear  flows  are  to  be  better  described,  entirely  new  models  are  needed 
for  the  deviatoric  part  of  the  pressure-strain  correlation  that  incorporate  some  compressible 
turbulence  physics. 
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Tkble  1.  Comparison  of  the  model  piedictions  for  the  equilibrium  Reynolds  stress  anisotropies 
with  the  DNS  results  of  Rogers  et  al  [23]  for  incompressible  homogeneous  shear  flow. 


Equilibrium  LRR  Model  FLT  Model  SSG  Model  DNS  Data 
Values 


-0.130 


-0.036 


0.189 


II 


-0.148  -0.165  -0.118 


-0.138  -0.148  -0.236 


-0.051  -0.082  -0.188 


'iDtble  2.  Comparison  of  the  compresribk  model  predictions  (using  the  dilatational  terms  of 
Sarkar  et  at  [14,  15])  for  the  equilibrium  Reyncdds  stress  anisotropies  with  the  DNS  results 
of  Blaisddl  et  ai  [3]  tot  the  compressible  homogeneous  shear  flow. 


Equilibrium 
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-0.137 

-0.148 

-0.236 

633 

-0.036 

-0.050 

-0.083 

-0.188 

IVible  3.  Comparison  of  the  compressible  model  predictions  (using  the  dilatational  terms 
of  Zeman  [16, 17])  for  the  equilibrium  Reynolds  stress  anisotropies  with  the  DNS  results  of 
Blaisdell  et  al  [3]  for  compressible  homogeneous  shear  flow. 
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Figure  1.  Time  evolution  of  the  turbulent  kinetic  energy:  Comparison  of  the  model  pre¬ 
dictions  (without  explicit  dilatational  terms)  and  the  DNS  results  of  Blaisdell  et  al  [3]  for 

compressible  homogeneous  shear  flow.  ( - )  LRR  Modd;  ( - )  FLT  Model;  ( - )  SSG 

Model;  (O)  DNS  results. 
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Figure  2.  Time  evolution  of  the  turbulent  diasipntion  rate:  Comparison  of  the  model  pre¬ 
dictions  (without  explicit  dilatational  terms)  and  the  DNS  results  of  Blaisdell  tt  al  [3]  for 

compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  ( - )  FLT  Model;  ( - )  SSG 

Model;  (O)  DNS  results. 
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figure  3.  Hme  evolution  of  the  turbulent  Much  number:  Comparison  of  the  model  pre¬ 
dictions  (without  explicit  dilatational  terms)  and  the  DNS  results  of  Blaisdell  et  al  [3]  for 

compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  ( - )  FLT  Model;  ( - )  SSG 

Model;  (O)  DNS  results. 
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ilgaxe  4.  Time  evolati<m  of  the  turbulent  kinetic  energy;  Comparison  of  the  model  predic¬ 
tions  (with  the  dilntntional  terms  of  Sarkar  e<  ol  [14,  15])  and  the  DNS  results  of  Blaisdell 

ei  al  [3]  for  compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  (•--•)  FLT  Model; 

( - )  SSG  Model;  (O)  DNS  results. 
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Figure  5.  Time  evolution  of  the  turbulent  duudpution  rate:  Comparuon  of  the  model  predic¬ 
tions  (with  the  dilatational  terms  of  Sarksi  et  ai  [14,  15])  and  the  DNS  results  of  Blaisdell 

et  ai  [3]  for  compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  (----)  FLT  Model; 

( - )  SSG  Model;  (O)  DNS  results. 
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Fignxe  6.  'nme  eWntion  of  the  tvibnlent  Mach  number:  Comparison  of  the  model  predic- 
ti<ms  (with  the  dOatatkmal  terms  of  Sarkar  e<  of  [14,  15])  and  the  DNS  results  of  Blaisdell 

e<  of  [3]  for  compresaihle  homofeneous  shear  flow.  ( - )  LRR  Model;  (--■-)  FLT  Model; 

( - )  SS6  Model;  (O)  DNS  results. 
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F^gnxe  7.  1^6  evolution  of  the  turbulent  Idnetic  energy:  Comparison  of  the  model  predic¬ 
tions  (with  the  dilatational  terms  of  Zeman  [16, 17])  and  the  DNS  results  of  Blaisdell  et  al  [3] 

for  compresdble  homogeneous  shear  flow.  ( - )  LRR  Model;  (----)  FLT  Model;  ( - -) 

SSG  Model;  (O)  DNS  results. 


Fignie  8.  Time  evolution  of  the  turbulent  dissipation  rate:  Comparison  of  the  model  pre¬ 
dictions  (with  the  dilatational  terms  of  Zeman  [16, 17])  and  the  DNS  results  of  Blaisdell  et 

ai  [3]  for  compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  (----)  FLT  Model; 

( - )  SSG  Model;  (O)  DNS  results. 
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Figure  9.  Time  evolution  of  the  turbulent  Mach  number:  Comparison  of  the  model  predic¬ 
tions  (with  the  dilatational  terms  of  Zeman  [16, 17])  and  the  DNS  results  of  BlaisdeU  et  ai  [3] 

for  compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  (----)  FLT  Model;  ( - ) 

SSG  Model;  (O)  DNS  results. 
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Figure  10.  Time  evolution  of  the  Reynolds  stress  anisotropies:  Comparison  of  the  model 
predictions  (with  the  dilatational  terms  of  Sarkar  et  al  [14,  15])  and  the  DNS  results  of 

Blaisdell  et  at  [3]  for  compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  (----) 

FLT  Model;  ( - )  SSG  Model;  (O)  DNS  results,  (a)  bji,  (b)  622  Mid  (c)  612* 
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Figure  10.  'Bme  evolution  of  the  Reynolds  stress  anisotropies:  Comparison  of  the  model 
predictions  (with  the  chlatational  terms  of  Sarkar  et  ai  [14,  15])  and  the  DNS  results  of 

Blaisdell  et  al  [3]  for  compressible  homogeneous  shear  flow.  ( - )  LRR  Model;  (----) 

FIT  Model;  ( - )  SSG  Model;  (O)  DNS  results,  (a)  5ii,  (b)  622  and  (c)  613. 
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